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Abstract 

Using the Robertson projection operator formalism, we derive generalized 
Bloch equations which describe the dynamics of a biased two-level tunneling 
system strongly driven by an external field and weakly coupled to a super- 



Ohmic heat bath. The generalized Bloch equations constitute a set of coupled 
nonlinear integro-differential equations. With their help we investigate the 
influence of phonons on the phenomenon of dynamical localization. 
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I. INTRODUCTION 



The last decade has witnessed detailed investigations of the problem of macroscopic 
quantum coherence .0'! The generic model is a two-level tunneling system (TLS) coupled to 
a harmonic heat bath and an external driving field. The tunneling units may represent, e. 
g., small groups of atoms or molecules in structurally disordered solids, such as in the study 
of low-temperature properties of glasses!; electrons in semiconductor quantum wells! or in 
the electron-transfer dynamics in chemical reactions^; or magnetic flux quanta at Josephson- 
junctions or SQUIDsi/ The basic mechanism by which quantum coherence is destroyed is the 
coupling of the quantum system to the dissipative environment. Various approaches using 
second-order perturbation theory, variational polaron theory^, mode-coupling theories^ or, 
most popular, the noninteracting-blip approximation (NIBA)Bi have been used to discuss 
the decoherence effect of the bath. 

One key issue is whether by careful tuning of the driving field the effect of dissipation 
may be reduced and quantum coherence be restored on the macroscopic level. This idea of 
controlling tunneling and relaxation has become especially popular since the discovery, by 
Grossmann, Hanggi and coworkers^, of the effect of coherent destruction of tunneling. In 
a numerical treatment of an isolated bistable quartic potential driven by a periodic force - 
but not coupled to a heat bath - they found complete suppression of tunneling: at certain 
ratios of the field amplitude and the driving frequency the two lowest Floquet eigenstates 
(quasienergies) of the driven system become degenerate, thereby preventing the system 
from performing coherent oscillations between the left and right well. As a result, a particle 
initially localized in one well will remain there forever. This so-called dynamical localization 
effect has also been described in a two-state modelJlUll and has since been the subject of 
several further investigationslliHil 

Recent studies have focused on the question to what extent the dynamical localization 
of a TLS will be affected by the coupling to a heat bath. Physical intuition seems to suggest 
that in order for dynamical localization to occur, a wave paket initially localized in one well 
must have a well-defined phase; only then can one expect that a mismatch between the 
tunneling motion and the driving field may prevent a particle from escaping into the other 
well. Since phonons destroy this phase coherence i. e., there is a finite phase memory time 
T2 - dynamical localization should always be softened through the coupling to a heat bath. 

Recently Grifoni et a/.lll gave a systematic approach to the transient and steady-state 
dynamics of a driven dissipative TLS. They considered Ohmic and frequency-dependent 
damping in the framework of the NIBA for the stochastic force, but without any approx- 
imation for the driving force. In Ref. [14| they extended their approach beyond the NIBA. 
Dittrich et investigated numerically the effect of driving and dissipation on the coherent 
tunneling motion of a symmetric bistable system for weak Ohmic damping. DakhnovskiSO 
employed small-polaron theory to address the same issue in a two-state approximation within 
the NIBA. He concluded for the case of a symmetric TLS with super-Ohmic damping that, 
in contrast to intuition, the localization transition remains stable against disturbance by the 
bath as long as the driving frequency is larger than the relaxation energy (polaron-binding 
energy) of the lattice. On these grounds he predicted the existence of a slow mode oscillation 
near the transitionJll His conclusions essentially rely on the fact that near the localization 
transition the rate of phase loss - which is due to the bath - decreases faster than the 
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tunneling coherence. 

So far most investigations have made use of the NIBA. Yet the NIBA breaks down at 
non-zero bias, weak dissipation and low temperatures. It is the latter regime, so far poorly 
understood, which we shall discuss in this paper. Weak TLS-phonon coupling is the limit 
which is relevant for, e. g., the description of tunneling centers in dielectric solids. We 
will work within the two-state picture and assume a super-Ohmic spectral density for the 
phonons; and we will treat the TLS-phonon coupling in Born approximation.il (There are 
in fact general arguments why under the above assumptions it should always be justified to 
treat the dissipation perturbatively.0) The external field, on the other hand, will be treated 
to any order. We are then able to derive nonlinear and non-Markovian generalized Block 
equations (GBE). These in turn permit us to discuss the influence of weak dissipation on the 
localization transition for both the symmetric and the biased case, and for arbitrarily strong 
driving and any temperature. In the high-temperature and high-frequency limit there are 
analytical solutions, which we will compare with those obtained from the NIBA. 

The paper is organized as follows. In Sec. II we give a short introduction to the Robertson 
formalism and apply it to driven TLS-dynamics. In Sec. Ill we derive the GBE. The linear 
response regime is discussed as a special case in Sec. IV. In Sec. V we address the influence 
of weak dissipation on the dynamical localization transition, before we close (in Sec. VI) 
with a brief summary. 



II. ROBERTSON FORMALISM 

We will derive a non-Markovian generalization of the Bloch equations with the help of the 
so-called Robertson formalism^. The Robertson formalism constitutes a particular variant 
of the well-known projection techniquei!, one that is tailored to describing the evolution of 
selected observables even far from equilibrium. 

When studying the dynamics of a macroscopic quantum system, one is typically con- 
fronted with the problem of determining the evolution of only a small set of selected ( "rel- 
evant") observables {G a }. Let the (generally time-dependent) Hamiltonian be denoted by 
H(t). Then the equation of motion for the expectation values 

g a (t) = (G a ) p(t y.= tT[p(t)G a ] (2.1) 

reads 

g a (t) = i(£(t)G a ) p{t) , (2.2) 

with p(t) being the statistical operator and C(t) the Liouvillian 

Cft^h-^Hit),*] . (2.3) 

Associated with the expectation values {g a (t)} is, at each time t, a generalized canonical 
state 

p Iel (t) := Zity 1 exp (- l a (t)G a ) , (2.4) 



2 



called the relevant part of the statistical operator, with partition function 



Z(t) :=tr exp f-£* a (t)GU 



(2.5) 



and the Lagrange parameters {l a (t)} adjusted such as to yield the correct expectation values 
( [2.1|) of the relevant observables. The difference 

Pirr(t) :=p(t) -Prel(t) (2.6) 

is then the irrelevant part of the state. 

Like all variants of the projection technique, the Robertson formalism is based on a clever 
insertion of projection operators into the equation of motion (|2.2|). Here the projection 



operator V(t) is chosen such that it projects arbitrary vectors in Liouville space onto the 
subspace spanned by the unit operator and by the relevant observables {G a }, and that this 
projection is orthogonal with respect to the (time-dependent) scalar product 



(A;B)®:= / d/itr M^W*) 1- ^ . (2.7) 



o 

Since this scalar product varies in time, the projection operator, too, is time-dependent; so is 
its complement Q(t) := 1 — V(t). The projection operator is known as the Kawasaki- Gunton 
projector^ 

Let T(t', t) be the evolution operator defined by 

| 7 T(t', t) = -i Q(t')C(t') Q(t')T(t', t) (2.8) 

with initial condition T(t,t) — 1. As Q projects out the irrelevant component of an observ- 
able, this operator may be pictured as describing the evolution of the irrelevant degrees of 
freedom of the system. The equation of motion for the expectation values of the relevant 
observables can then be cast into the form 

g a (t)=i(C(t)G a ) iem + f &t'Y.{Q(t ! )C{t')G c -T{t\t)Q(t)C{t)G a fn c (t') 

Jo c 

+i(T(0,t)Q(t)£(t)G a ) irm , (2.9) 

which is known as the Robertson equation^ Here (-) re i and (■){„ denote expectation values 
evaluated in the relevant and irrelevant part of the state, respectively. A major simplification 
occurs if, as is often the case, the initial macrostate can be characterized completely by 
the expectation values {(G a )(0)} of the relevant observables and hence has the generalized 



canonical form ( |2.4| ). In this case the initial state has no irrelevant component, which in 
turn implies that the third ( "residual force" ) term vanishes. 

The Robertson equation embodies a closed system of nonlinear coupled integro- 
differential equations for the expectation values {g a (t)}. These coupled equations are non- 
local in time: future expectation values of the relevant observables are predicted not just on 
the basis of their present values, but on their entire history. In essence, the irrelevant degrees 
of freedom have been eliminated, but this elimination must be paid for by nonlinear and 
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non-Markovian features of the projected equation of motion; the complexity of the evolution 
has been mapped onto nonlinear it ies and a non-local behavior in time. 

We now determine the Robertson equation for a driven TLS. The underlying microscopic 
dynamics of the TLS is determined by a time-dependent Hamiltonian 

H(t) = H s (t) + H SB + H B . (2.10) 

Here the system Hamiltonian is given by 

H s (t) = -\a-h(t) , (2.11) 

where h(t) = % ■ (A , 0, —A — /(£)) includes the tunneling frequency A , the bias A, and the 
time-dependent driving field 

f(t) = 2ficos(cj L t) (2.12) 

with Rabi frequency Q and driving frequency loj_,. If we denote the state of a particle localized 
in the left (right) well by \L) (\R)) then the Pauli matrix a z = \L)(L\ — \R)(R\ represents its 
coordinate, while a x — \L) (R\ + \R)(L\ represents its tunneling motion. With this convention 
A > means that the left well has a higher potential energy than the right well. The bath 
Hamiltonian H B = Y,k frwkb\bk, with bosonic creation and annihilation operators satisfying 
[bi, bj] = Sij, describes the dynamics of an ensemble of harmonic oscillators (phonons). Their 
spectral density is taken to be super-Ohmic, 

2 

J(co) = -y c?5(u — ujj) = Uco 3 exp(— uj/u c ) , uj > , (2-13) 

up to some cutoff u c . Finally, the system-bath coupling is assumed to be bilinear in the 
tunneling and the phonon coordinate, 

H SB = a z J2c J (b J +b})=:a z e . (2.14) 

j 

Relevant observables are 



{G a }^{a a ,(H SB + H B )} , (2.15) 

with expectation values 

{g a (t)}^{ Pa (t),E B (t)} (2.16) 

and associated Lagrange parameters 

{I a (t)}^{-j3\a(t)/2,j3} . (2.17) 

Here p a (t) denotes the time-dependent polarization of the TLS, and E B (t) the internal energy 
of the bath. We assume the heat bath to be large enough so that its inverse temperature 
(3 = 1/ksT does not vary in time. The relevant part of the statistical operator reads 
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Pvel{t) 



Z exp 
Z~ x exp 



-P (H B + H SB - ±A(fj • B 
-(3 (H(t) + \${t)-\{t))-B 



(2.18) 



Provided the (full) initial state p(0) can be written in this generalized canonical form, i. e., 
provided p(0) = p re i(0), then the residual force term in the Robertson equation vanishes. 
The class of initial states in the generalized canonical form includes the common case of a 
particle initially held at the site \L), p z (0) = 1, and coupled to a bath in thermal equilibrium, 



p(0) oc 7T L exp [-(3 (H(0) + l%A a 3 



7T L 



(2.19) 



Here, 7Tl = |(1 + er 2 ) denotes the projector on the state \L). This initial condition has 



the generalized canonical form ( |2.4| ) with Lagrange parameters A(0) = (0,0, oo). Clearly, 
the fact that the initial state has the generalized canonical form does not imply that the 
total system is in equilibrium; rather, it means that the initial state can be completely 
characterized by the expectation values {p a (0), E B (0)} of the relevant observables. The 
Robertson equation then translates into the generalized Bloch equations (GBE) 



Pa(t) =n- L [p(t) X (h(t) - X(t))] a 



P 



with memory kernel 



K aP (t,t') := §(Q(t')C SB (Tp;T(t',t)Q(t)C SB a a ) 



(*') 



(2.20) 



(2.21) 



The GBE are non-Markovian: our formalism allows for the inclusion of arbitrarily large 
memory effects. Indeed, such non-Markovian effects will prove essential for the determination 
of the correct relaxation rates. 



III. DRIVEN DYNAMICS AT WEAK DISSIPATION 

For simplicity, we will from now on set % = ks = 1- In order to evaluate the memory 
kernel we assume weak coupling between the TLS and the bath, and hence do lowest (i. e., 
second) order perturbation theory in H$ B (Born approximation). To this order we can omit 
the complement projectors Q, replace the evolution operator T with 

T(t',t) -> U s (t',t)®U B (t',t) , (3.1) 

where T, Us, U B are the evolution operators associated with QCQ, and C B , respectively, 
and evaluate the scalar product (•; •p t ' in the state 

pSW =PB®Ps[A(t)] 

:= exp(-(3H B ) ® ±- exp(|/?A(t) • B) (3.2) 

rather than in p Te \(t). We thus obtain 
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K af3 (t,t')=g(t-t')-E a(S (t,t') 



(3.3) 



with the bath relaxation function 

g(t - t') := 2(3(e; U B (t' , t)e) B = 2 H doj cos(uj(t - t')) (3.4) 

Jo UJ 

and spin relaxation function 

S Q/3 (M') :=\([a z ,a p ]-U s (t\t)[a z ,a a ]) { p . (3.5) 

Here (•;•)!* and (•; -) B denote scalar products evaluated in the states p s [X(t')] and p B , 
respectively. As the spectral density J(ou) has a characteristic width oj c , the bath relaxation 
function g{r) -which is essentially its Fourier transform- decays on a typical scale ("memory 
time") r mcm ~ 1/uJc- For super-Ohmic damping this memory time is the only characteristic 
time scale of the bath, independent of the temperature. This distinguishes the super-Ohmic 
from the Ohmic case where there exists a second characteristic scale of order (KT)^ 1 , where 
K is the Kondo parameter. 

Using now, in Born approximation, 

m = xit'W) (3-6) 

with 

x(0-^tanh^M , (3 .7) 
V ; |A(t')| 2 

and defining the instantaneous equilibrium polarization 

(<?Ut>) ■= X(t')h(t') , (3-8) 
the GBE can be cast into the compact form 

p a (t) = \p(t) x h(t)] a - f dt'g{t-t')Y,^ap{tJ)x~\t')\pp{t')-(<yp)h{t')\ ■ (3.9) 

This approximate equation of motion, valid at weak dissipation, is still non-Markovian, and 
still nonlinear in p. Far from equilibrium the nonlinearities may become significant and lead 
to a non-exponential relaxation of the polarization vector. Also the external driving field 
enters nonlinearly, through the time evolution operator in the spin relaxation function. 



IV. LINEAR RESPONSE REGIME 

To test our generalized Bloch equations we show that in the linear response regime they 
yield results consistent with earlier calculations. We assume that the dynamics take place 
in the linear regime, i. e., that at all times the polarization be close to its instantaneous 
equilibrium, p ^ (a). In this regime we may evaluate both H(t,t') and m the 
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instantaneous equilibrium state p s [h(t')], rather than in the state ps[X(t')}, thus rendering 
the GBE linear in p. Approximating the spin relaxation function by 



/ cos(e r) — Uosin(e r) 
M sin(e r) ul + Mq cos(e r) | (4.1) 

V o oo 

with r :—t — t', eo := \J A 2 + Aq, u := A /e and u := A/e , and expanding the suscepti- 
bility in a Fourier series 

oo 

X-\t)= £ x^e" TOi , (4.2) 

m=— oo 

the GBE can be solved by Laplace transformation. Keeping terms up to linear order in Q 
one finds, after tedious but straightforward calculations, 

PM=pfit)+p^\t) (4.3) 

where the first term describes the transient dynamics, while the second term describes the 
steady-state dynamics. The transient term is given by the standard second order result for 
a static bias,E3 

rf 0) (*) = 4 e- t/T 2 0) cos(e t - (p)l cos <p + (u 2 - p^) + p^ , (4.4) 

with p e z q = — u tanh(/3e /2) and coty? = e^r^, and relaxation rates 



IT. 



(0)\-l _ /o,J0)\-l _ 1„ 2 77J 



(2rr)- 1 = X7r^coth(/3e /2) . (4.5) 



The steady-state term gives the linear response of the TLS to the driving field. In the 
asymptotic regime where all transients have decayed, the polarization oscillates coherently 
with the driving frequency Ul; in the linear regime there is no generation of higher harmonics. 
Defining the dynamical susceptibility x( u l) via the relation 

= -2n (xM * + x(-uJl) e Wht ) (4.6) 
we find for resonant driving (uj l « ±e ) 

i- (o) 

kJul) = | tanh(/3 eo /2) £ : ., ~ - m (4-7) 

± 1 - i{ujl =F eo) r 2 v 



and for low-frequency driving (cj^ <C e ) 



0/4 w 



n2 



Xrci(^L) - Xoo = , 2 ,„ , ; , (4.8) 

cosh (/?e /2) l _ %UlT \ u > 

where Xoo = ( M o/4 e o) tanh(/3e /2) denotes the instantaneous response. In the derivation we 
have made use of XoX±i — ~ {P e o coth(/3e /2))/(2 cosh 2 (/3e /2))]. These results are 

well-known from the literature, see e. g. Ref. 
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V. EFFECT OF WEAK DISSIPATION ON DYNAMICAL LOCALIZATION 



We return to the general case of dynamics arbitrarily far from equilibrium. We factor out 
that part of the dynamics which is due to the driving field alone, by defining new polarization 
vectors 

P 



(? Q >ft(i):=E^(i)(^)Mf) (5.2) 



via an orthogonal rotation matrix 

R(t) := 



/ cosF(t) sinF(t) 

-sinF(t) cosF(t) | (5.3) 
V 1 



with 

F(t):= / f(s)ds = —sm(u L t) . (5.4) 
Jo uo L 

This rotation affects only the x- and ^/-components of pit), but not p z (t). The thus rotated 
polarization vectors obey an equation of motion in which the driving field no longer appears 
explicitly, but only indirectly via R(t); the only field which still appears explicitly is the 
static part 



h st ,f.= h(t) =(A ,0,-A) . (5.5) 
The spin relaxation function in the GBE must be rotated, too, and is now given by 

(i2(t)S(t,Oi2 _1 (0)a/9= ^[^^jWn^OWs^O^^O)^,^])^ , (5.6) 

where U Q denotes the spin evolution operator associated with the driving field alone, at zero 
static field. Since 

(R(t)E(t,t')R-\t r )) zl = \Ji = x,y,z , (5.7) 
the ^-component of the polarization vector satisfies the equation of motion 

jp z (t) = -A [cos F{t)p y {t) + sin F{t)p x {t)\ . (5.8) 

This equation, together with the corresponding (more complicated) equations for p x (t) 
and p y (t), has provided the basis for an earlier analysis^ of the unbiased (A = 0), dissipa- 
tionless TLS. In Ref. [12| the authors find in the fast-driving limit (u>l 3> Ao) the driven TLS 
to be equivalent to a TLS without external field, but with renormalized tunneling parameter 
A J (2Q/u L ) instead of A ; and, accordingly, describe the dynamics by 
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p z (t) = cos[A J (2fi/cu L )t] 



(5.9) 



The authors conclude that in the fast-driving limit and for small but nonzero Jq{2VL/ujl) 
the dynamics exhibits strong slow mode oscillations ( "low- frequency generation"); while for 
J (2Q/ujl) = the component p z (t) stops evolving at all: it remains forever "frozen" in its 
initial state ("dynamical localization"). 

In order to understand and later generalize this result in our formalism we, too, invoke 
the fast-driving limit: we assume that the driving frequency be much larger than both 
tunneling frequency and bias, 

w L >A ,A. (5.10) 

Yet we do not assume the TLS to be unbiased or dissipationless. The fast-driving limit 
allows one to replace rapidly oscillating terms by their time averages: 



R(t)h stat -> R(t)h stat = {A J {2Q/u L ), 0, -A) =: h eS (5.11) 

@)h(?) -> x(t')R(t')h(t') =: p as (5.12) 

(R(t)E(t,t')R-\t')U - Wm^tWW))t(s =■ Z a p(t - • (5.13) 

The latter average is taken at fixed (t — t r ). With these replacements the GBE (|3.9|) acquires 
the convolution form 

P*(t) = $(t) x h eS ] a - f dt'g(t-t')Y,^(t-t')x-\t ! )ipp(t')- P ™} ■ (5.14) 
in particular, 

jMt) = -J (2Q/u L ) A p y (t) (5.15) 

which predicts a slowing down of the time evolution of p z (t) near Jo = 0. 

To this time-averaged GBE there are rapidly oscillating correction terms which are neg- 
ligible as long as Jo 7^ 0, but which may become important at the localization transition 
Jb = 0. Infed, numerical studies by Makarov and MakriS indicate that at Jo = the 
rapidly oscillating correction terms cause additional phase relaxation (if the system is coupled 
to a heat bath), and eventually destroy the dynamical localization. Hence the time-averaged 
GBE ( |5.14| ) can only describe the dynamics near, but not exactly at, the localization tran- 
sition. Furthermore, even away from the localization transition, the time-averaged GBE 
misses short time transients (Gaussian or algebraic decay) arising from the bath decorrela- 
tion and from the nonlinear time evolution induced by the laser field. Our description is 
thus expected to be valid on intermediate time scales uj^.ujJ} <C t r a , where the {t^ 1 } 
are the rates at which the various components of the polarization vector p(t) relax to their 
stationary state. 

A full analytical solution of the dynamics of the polarization vector p(t), and hence a 
study of its dynamics close to the roots of the Bessel function J , is only possible if we 
invoke further approximations: namely, the high-temperature limit, T ^> A , A,f2, and the 
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high-frequency limit, ujl ^> u c . In the high-temperature limit we may evaluate both E(t, t') 
and i n Ps[0] = |ls rather than in the time-dependent state p s [\(t')}, allowing us to 

replace 

(•;->r - (-;-}o , x - 1 (t')^2/p . (5.16) 



The product of evolution operators in the rotated spin relaxation function ( |5.6| ) can be 
written as a time-ordered exponential 



U n (0,t')U s (t',t)U Q (t,0) = Texp 



% 

2 Jt> 



fdt"{R{t")h stat )-CT* 

Jt> 



(5.17) 



where cr x denotes the commutator with a. In the high-frequency limit the integrand in 
the exponent oscillates rapidly during the interval [t',t], which has a typical length r mom ~ 
1/oj c ^> 1/oJl. In order to perform the average ( [5 . 1 3| ) we may thus replace 

f dt"(R(t")h stat ) — > (t - f) ■ R(t")h s J , (5.18) 

which will yield just the time evolution operator of a static tunneling problem with new 
parameters (|5.11|) . Defining 



e:=^A 2 + A 2 J 2 (2Q/uj L ) , (5.19) 
u := A J (2fi/w L )/e , (5.20) 
u ■= A/e = Vl - u 2 , (5.21) 



the rotated spin relaxation function is then given by 



/ cos(er) — -usinfer) 



a t 



wsin er u 



2 



w 2 cos(er) . (5.22) 



V 



We recognize that now the dynamics has the same form as in the linear response regime 
(Sec. IV), with the additional simplification /?e 1. Hence, indeed, the driven system can 
be mapped onto a time-independent TLS with modified parameters. The results for the 
time evolution of p z (t) can immediately be taken over from (|4.4j ), with the sole replacement 

A -> A J (2n/uj L ) ; (5.23) 

and consequently, u — > m, u — > u and e — > e. Noting that (t^)- 1 oc A 2 (Eq. |3) we find 
that the relaxation rates must be modified according to 

(t^^IM^/^)] 2 ^)- 1 ■ (5-24) 

In the unbiased case (A = 0) and away from the roots of the Bessel function the system 
thus reaches its steady state on the time scale /[J (2Q/ujl)} 2 . Near the localization 
transition (J — * 0) this time scale diverges faster than the tunneling time (A Jq{2VL/ujl))~ 1 , 
dissipation thus being effectively switched off. As one reads off from the linear response 
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solution ( |4.4| ), this implies that there will be a slow mode coherent oscillation near the 
localization transition, and that the particle will remain trapped in one well on the time 
scale (A Jq(2Q/uji j ))^ 1 . This is consistent with Dakhnovskii's earlier analysis within the 
framework of the nibaB and with the ab initio numerical calculations of Makarov and 
MakrS 

In the biased case (A ^ 0), on the other hand, the amplitude of the oscillatory term in 

(ED, 

u 2 = (J (2fi/^)A ) 2 /e 2 , (5.25) 

becomes negligible for J (2Q/uj l ) <C A/A . One thus expects pure exponential decay, rather 
than slow mode oscillations, in the immediate vicinity of Jo = 0. (One way of viewing this 
phenomenon is that as the bias is switched on, the quasienergy levels cease to cross and 
hence the destruction of tunneling can no longer be coherent.) The particle will then remain 
trapped in one well on the time scale r[ /[Jo(20/o;i)] 2 . Further away from the localization 
transition the exponential decay will be superimposed with small-amplitude oscillations. 

The evolution of p z (t) in the two cases (unbiased and biased) is illustrated in Figs. 1 and 
2, respectively. We have calculated the curves using Eq. (|4.4j) , with the replacement ( |5.23|) . 



VI. SUMMARY AND FINAL REMARKS 

In this paper we have studied the influence of a weakly coupled thermal environment 
on dynamical localization. Motivated by the breakdown of the NIBA at non-zero bias and 
weak dissipation we have investigated the dynamics of a biased two-level tunneling system 
strongly driven by an external field and weakly coupled to a super-Ohmic heat bath. We 
have employed the Robertson projection operator formalism to derive generalized Bloch 
equations for the polarization vector p(t). 

Assuming only the validity of the two-state picture and of the Born approximation, and 
invoking the fast-driving, high-frequency and high-temperature limits, we have shown that 
on intermediate time scales the driven dissipative TLS can be mapped onto a dissipative 
TLS without driving, but with renormalized parameters. We have found that since the 
phase relaxation rate decreases quadratically with Jq{2VL/uil), while the tunneling frequency 
decreases only linearly, coherence is restored near the localization point of an unbiased TLS 
(A = 0). This manifests itself in slow mode coherent oscillations near the localization 
transition, a result that Dakhnovskii had previously obtained within the framework of the 
NIBA. In the biased case (A ^ 0) we found a qualitatively different behavior, namely pure 
exponential decay in the immediate vicinity of Jo(2Q/ujl) = 0. 

Let us finally comment on two other papers which address the problem of driven tunnel- 
ing dynamics in the weak-dissipation limit 00 Dittrich, Hanggi and OelschlagerEl consider 
dynamical localization in the weak-coupling limit (Born approximation) and solve numer- 
ically the quantum master equation for the full continuous bistable system with Ohmic 
dissipation. In contrast to our present analysis, however, they employ a restricted rotating- 
wave approximation and assume the bath to be Markovian, i. e., they assume that the 
correlation time for the boson modes is negligible compared to the characteristic time scale 
of the double-well dynamics. For Ohmic friction this assumption is justified if in addition to 
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the classical time scale r mem ~ l/ou c also the quantum time scale of the bath (KT)" 1 , with 
K being the Kondo parameter, is shorter than the tunneling time. But in this case the dy- 
namics is incoherent and a weak-coupling picture inappropriate. In fact, it has been shown 
explicitly that the Born approximation fails for a tunneling system subject to Ohmic dis- 
sipation and driving, and that driving renders the dynamics intrinsically non-Markovian.0 
This is not surprising since even for super-Ohmic dissipation, where the second (quantum) 
time scale does not exist, non-Markovian effects are essential for the determination of the 
correct relaxation rates.Ba 

Very recently Grifoni et a/.0 have extended their previous treatment^! beyond the NIBA. 
They derive an exact master equation, which they then use to determine the dynamics at 
weak coupling and fast driving. Indeed, our GBE finds a counterpart in the set of equations 
adjacent to their Eq. (10). Yet instead of discussing the influence of weak dissipation on 
the dynamics near the localization transition, they use their master equation as a starting 
point to investigate the modification of quantum coherence away from the roots of the Bessel 
functions, in the regime VL/ujl <C 1. 
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FIGURE CAPTIONS 

FIG. 1 Time evolution of p z (t) in the symmetric case. The parameters are T = 10A , Aqt^ - 1 = 
10 and A = 0, for driving fields with 2VL/uj l = ( — ), 2VL/uj l = 2 (- -), and 2VL/uj l = 
2.3 ( — ). The localization transition occurs at 2VL/uo L rs 2.405. 

FIG. 2 Time evolution of p z (t) in the biased case. The parameters are chosen as in Fig. 1, 
except for A = 0.75A . 
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